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ABSTRACT 

We present bright galaxy number counts in five broad bands (u', g' , r', i', z') from 
imaging data taken during the commissioning phase of the Sloan Digital Sky Survey 
(SDSS). The counts are derived from two independent stripes of imaging scans along 
the Celestial Equator, one each toward the North and the South Galactic cap, covering 
about 230 and 210 square degrees, respectively. A careful study is made to verify the 
reliability of the photometric catalog. For galaxies brighter than r* = 16 , the catalog 
produced by automated software is examined against eye inspection of all objects. Sta- 
tistically meaningful results on the galaxy counts are obtained in the magnitude range 
12 < r* < 21, using a sample of 900,000 galaxies. The counts from the two stripes 
differ by about 30% at magnitudes brighter than r* = 15.5, consistent with a local 
2a fluctuation due to large scale structure in the galaxy distribution. The shape of 
the number counts-magnitude relation brighter than r* = 16 is well characterized by 
N oc 10°' 6m , the relation expected for a homogeneous galaxy distribution in a "Eu- 
clidean" universe. In the magnitude range 16 < r* < 21, the galaxy counts from both 
stripes agree very well, and follow the prediction of the no-evolution model, although 
the data do not exclude a small amount of evolution. We use empirically determined 
color transformations to derive the galaxy number counts in the B and Igu bands. We 
compute the luminosity density of the universe at zero redshift in the five SDSS bands 
and in the B band. We find £b = 2.4 ± 0.4 x lO 8 L /iMpc -3 , for a reasonably wide 
range of parameters of the Schechter luminosity function in the B band. 

Subject headings: cosmology: observations 



1. Introduction 

The number counts of galaxies as a function of magnitude is one of the classical cosmological 
tests. It has been repeatedly studied over many years by many authors, and has yielded valuable 
insights about the properties of galaxies and their evolution with redshift (see Sandage (1988) 
and Koo & Kron (1992) for reviews). However, despite much effort, the normalization of the 
counts-magnitude relation for bright galaxies remains uncertain by as much as 50%. This leads to 
uncertainties in the normalization of the galaxy luminosity function at zero redshift at the same 
level, and therefore to uncertainties in the interpretation of data on galaxy evolution with redshift. 
The slope of the galaxy counts at bright magnitudes is also a matter of debate; some authors have 
claimed that it is steeper than the slope predicted by the no-evolution model (Maddox et al. (1990), 



1 Based on observations obtained with the Sloan Digital Sky Survey. 
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but see also Weir et al. (1995)), and argued either for a local underdensity in the galaxy distribution 
or for a rapid evolution of the galaxy population since z ~ 0.2. 

Bright galaxies are rare on the sky, and their number density varies due to large-scale structure. 
Hence, we require imaging surveys extending over large solid angles to obtain reliable galaxy counts 
at the bright end. For this reason, most bright galaxy counts have been based on photographic 
plate material from Schmidt telescopes (Maddox et al. 1990; Bertin & Dennefeld 1997), for which 
it is notoriously difficult to obtain accurate photometric calibration. Surveys with CCD imagers 
(Hall & Mackay 1984; Tyson 1988; Lilly, Cowie, & Gardner 1991; Metcalfe et al. 1991) have made 
dramatic progress in determining the galaxy counts at faint magnitudes, but the field of view 
covered by these surveys is usually small, making it difficult to accurately measure the number 
counts at bright magnitudes. Even the ambitious surveys of Gardner et al. (1996) and Postman et 
al. (1998) are limited to 8.5 deg 2 and 14.7 deg 2 , respectively, much smaller than the fields surveyed 
using photographic plates. What is clearly needed is an imaging survey of the sky using modern 
CCD detectors, covering at least several hundred square degrees. 

The Sloan Digital Sky Survey (SDSS; York et al. (2000)) consists of an imaging survey in five 
photometric bands of -it steradians of the Northern sky, as well as a follow-up spectroscopic survey 
of roughly 10 6 galaxies and 10 5 quasars, complete within precisely defined selection criteria. Images 
of the sky are obtained by SDSS at the rate of 20 square degrees an hour. Thus, useful data to 
determine the galaxy counts at bright magnitudes may be obtained in a few nights of imaging. 
In this paper, we present galaxy number counts in the five SDSS passbands u', g 1 , r', i', and z', 
using imaging scans taken during the commissioning phase of the survey. We derive the galaxy 
counts in the magnitude range 12.5 < r* < 21, from two patches of the sky along the Celestial 
Equator, one each toward the North and the South Galactic cap, covering an area of about 230 
deg 2 and 210 deg 2 , respectively. We emphasize that the number counts of galaxies constitute a 
crucial verification of the photometric catalog derived from the imaging survey, in particular on the 
uniformity of photometric measurements with magnitude. 

The outline of this paper is as follows. In §2, we give a brief description of the observations and 
the data. We describe the photometry of galaxies, and the associated errors, in §3. In §4, we describe 
the construction of the galaxy catalog in the Northern equatorial stripe, from the full photometric 
catalog. We present the galaxy number counts in the five SDSS bands in §5. In §6, we convert 
our number counts to the B and I^ia bands using empirical color transformations, and compare 
them with existing galaxy count results in the literature. We calculate the luminosity density of 
the universe in the SDSS bands and in the B band in §7. As a consistency test we recalculate 
the normalization of the luminosity functions for the SDSS passbands, which are obtained from 
the spectroscopic survey of the SDSS project, and are published prior to this paper (Blanton et 
al. 2000). In §8, we derive the galaxy counts from the Southern equatorial stripe. We present the 
color distributions of galaxies in §9. We present our conclusions in §10. 
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2. Imaging Observations and Data 

The SDSS survey is carried out using a wide-field 2.5 m telescope, a large format imaging 
camera (Gunn et al. 1998), two fiber-fed double spectrographs, and an auxiliary 0.5 m telescope for 
photometric calibration. The sky is imaged in five passbands, u', g', r', i', z' covering the entire 
optical range from the atmospheric ultraviolet cutoff in the blue to the sensitivity limit of silicon 
CCDs in the red (Fukugita et al. 1996). The imaging camera contains thirty photometric CCDs, 
arranged in six columns of five rows, each row corresponding to a different filter. The data are taken 
in time-delay-and-integrate (TDI, or drift-scan) mode at the sidereal rate along great circles on the 
sky, yielding a strip consisting of six very long and narrow scanlines, each 13.5' wide. The effective 
exposure time (i.e., the transit time of an object across a single chip) is 54.1 sec. The scanlines of 
a given strip do not overlap, but observing a second strip offset from the first strip by about 12.8' 
gives a filled stripe 2.5° wide. The scanline in five colors produced by each column of five CCDs 
is divided into 2048 x 1489 pixel frames with an overlap of 128 pixels between adjacent frames; a 
set of five frames covering a given region of sky (in the colors u', g' , r', i' and z') is referred to as a 
field. The pixels are 0.396" square on a side, which satisfies the Nyquist limit for seeing > 1.0" full 
width at half maximum (FWHM) (see Gunn et al. (1998) for technical details). The 0.5m telescope 
observes standard stars to determine the atmospheric extinction on each night; the zeropoint of 
the 2.5m telescope is determined by observing patches of sky that the 2.5m is observing as well. 

The scans of the Northern equatorial stripe were obtained in two nights of SDSS commissioning 
observations on March 20 and 21, 1999, with the telescope parked on the Celestial Equator (SDSS 
Runs 752 and 756). The total observation time was about seven hours in Run 752 and about eight 
hours in Run 756. The two runs together form a filled stripe. For r' < 16, we present the galaxy 
counts from data in the right ascension range over which a full stripe is completed by this pair of 
runs: 145.°14 < a < 236.°42 and -l.°26 < 5 < l.°26 (all coordinates are in J2000), covering a 
total area of 230 deg 2 . 

The imaging data used in this paper were taken before the commissioning of the current 
Photometric Telescope. Therefore, we calibrated the data by observing secondary patches in the 
survey area using a (now decommissioned) 61 cm telescope at the observatory site, and by observ- 
ing primary standard stars using the US Naval Observatory's 40" telescope, with filter and CCD 
characteristics nominally identical to those at the SDSS Photometric Telescope. Since the transfor- 
mation from the primary standard stars to the objects observed with the SDSS 2.5m telescope has 
not yet been fully defined, we expect photometric errors of «0.05 mag in each band with respect 
to the proper SDSS photometric system. Thus, in this paper, we will denote our magnitudes as 
u*, g* , r*, i* , and z* , to emphasize the preliminary nature of our photometric calibration, rather 
than the notation u' , g' , r', i' , and z' that will be used for the final SDSS photometric system. 
However, we will use the latter notation to refer to the SDSS photometric passbands themselves. 
Our magnitude system is based on the AB95 system (Fukugita et al. 1996). The mean flux density 
over each broad passband is /„ = 3631 x lQ~ 0Am Jy. 
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Figure 1 shows the FWHM of the point spread function (PSF) in the r' band (the seeing) 
in each field of the third CCD column. The PSF is fit to a pair of concentric circular Gaussians; 
the FWHM shown in this figure is that of this model 23 . The top and bottom panels show the 
seeing variations in Runs 752 and 756, respectively. The median PSF FWHM in the r' band was 
1.4" in Run 752, and 1.3" in Run 756, although the seeing varied substantially over the duration 
of the runs. The sky brightness also varied significantly during the runs. For example in Run 756 
brightness of the v! and g' bands varied by 0.75 mag, that of r' by 0.45 mag, i' by 0.23 mag and z' 
by 0.61 mag during the scan for 8 hours. For the z' band we recognize a 10% variation in a short 
time scale (~ 10 min) in addition to a global change. The average sky brightness was 22.77, 22.06, 
21.05, 20.34, and 19.18 mag per square arcsec in u', g', r', i' , and z', respectively. 



3. Photometry of Galaxies 

The imaging data are processed with the photometric pipeline (hereafter Photo; Lupton et al. 
unpublished, see also York et al. (2000) and http://astro.princeton.edu/PBODK/welcome.htm) 
specifically written for reducing the SDSS data. Each field is processed independently by Photo. 
Galaxies do not have sharp edges, nor unique profiles of surface brightness distributions, which 
makes it non-trivial to define a flux for each object. Photo measures a number of different types 
of magnitudes for each object. We present the galaxy counts using the Petrosian flux (Petrosian 
1976), which is defined by 

/•krp 

F P = 2vr / I{r)rdr , (1) 



Jo 

where I(r) is the surface brightness profile of the object in question, and rp is the Petrosian radius 
satisfying 

„ = !M (2 ) 

' 27r/ rp I(r)rdr/(irr 2 P ) ' 1 } 

In practice, this implicit equation for rp is replaced with 

_ 2tt/ ^ I(r)rdr/ [vr((1.25r P ) 2 - (0.8r P ) 2 )] 
71 27r/ rp I(r)rdr/(nr%) ' 1 ' 

which yields robust measurements. Petrosian radius is independent of the distance of a galaxy and 
foreground reddening. Note that two parameters, k and rj, are needed to specify the Petrosian 
flux. We adopt k = 2 and rj = 0.2 to optimize requirements. A small r] is preferred in order to 
minimize the discrepancy in apertures for galaxies with the exponential profile and those with the 
de Vaucouleurs profile. A small r\ also makes the measurement of the Petrosian radius insensitive 
to seeing variations. On the other hand, consideration of signal to noise leads us to choose r\ > 0.2. 



23 Note that this FWHM measure gives a value somewhat smaller than that reported in some other SDSS papers, 
where the size of seeing is defined by the square root of the effective number of pixels calculated from the ratio of the 
first moment square to the second moment of the flux. 
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For every object, Photo calculates the Petrosian radius from aperture photometry using spline 
interpolation. With our adopted values of k and 77, the Petrosian flux is an integration over an 
aperture of radius of 3.6 half-light radii for an object with a de Vaucouleurs profile, and 7.2 scale 
lengths for an object with an exponential profile. To quantify the difference between Petrosian 
and total magnitudes for different galaxy morphologies, we carry out simulations with an empirical 
model for the mix of galaxy morphology and the scale size (Shimasaku et al. (2001)). The source 
of the other parameters used is explained in Section 5 below. We assume that the de Vaucouleurs 
profile is truncated at 5 half light radii 24 . We find that the Petrosian magnitude underestimates 
the total flux by ~ 0.03 ±0.01 mag (mean value) for this morphological mix, although the offset can 
be as large as 0.15—0.2 mag for face-on de Vaucouleurs profiles. The Petrosian radius is computed 
in all color bands, but we measure the Petrosian flux in any color band using the Petrosian radius 
for the r' band. 

We have compared the photometry of objects that lie in the overlap regions between Runs 752 
and 756. We find that the difference between the magnitudes of the same objects in these two runs 
taken under different seeing conditions is consistent with the errors in the magnitudes quoted by 
Photo. At m = 19, the errors in Petrosian magnitudes are about ±(0.03 — 0.04) mag in the g' , r', i' 
bands and about 0.05 mag in the u' and z' band. 

We have studied the completeness limit of the SDSS imaging data, by comparison with simu- 
lated images, and from comparison with deep HST images. Lupton et al. (unpublished) present a 
detailed study of the completeness limit of the SDSS imaging data. The comparison with deeper 
HST images shows that the SDSS imaging data is 50% complete to r* = 22.5, which agrees with the 
expectation from the simulation described above. This limiting magnitude is substantially fainter 
than the r* = 21 limit to which we present the galaxy counts in this paper. 

For the detection of objects Photo adopts a peak finding algorithm using the matched filter 
method, not only with the use of the PSF filter but also it is applied to 2 x 2 and 4x4 binned 
image. This algorithm is capable of detecting objects with considerably low surface brightness. Our 
simulation shows that the 100% detection completeness is sustained at r* = 21mag for galaxies 
with the scale length 5a larger than that of M* galaxies. 

4. Construction of the Galaxy catalog 

We have adopted two approaches to define a galaxy sample from the photometric catalog 
output by Photo. For bright magnitudes (r* < 16), we have visually classified all objects in Runs 
752 and 756 into stars and galaxies. In §4.3, we compare the results of this visual classification 



24 With the untruncated de Vaucouleurs profile, 11% of the flux is distributed outside the aperture of this radius. 
The total flux of simulations is renormalized to the truncated de Vaucouleurs profile. Without the truncation the 
tail of the de Vaucouleurs profile causes a 2—3% error in the local sky estimate. 
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with the results of the star-galaxy classification algorithm used in Photo. At fainter magnitudes 
(16 < r* < 21), we use the star-galaxy classification employed by Photo. In these faint magnitudes 
Photo yields consistent classifications of objects in different colors. 



4.1. Selection of objects from the photometric catalog 

Objects that lie in the overlap between adjacent scanlines in two strips of a stripe, and those 
that lie in the overlap between adjacent frames, appear more than once in the photometric catalog. 
Each of these detections is flagged as PRIMARY or SECONDARY, based on its position in the individual 
frames; the detections flagged as PRIMARY define a unique detection of each object. Therefore, we 
construct the galaxy catalog using only the PRIMARY objects. 

Photo examines every detected object for multiple peaks; if they exist, the object is deblended 
(Lupton et al. unpublished) into a hierarchy of component objects (children), and the parent is 
flagged as BLENDED. If a parent object has a large number of peaks, then only the 30 highest peaks 
are deblended. If the galaxy touches the edge of a frame, or if it is larger than half the size of a 
frame, it is not deblended any further, and the object is additionally flagged as NODEBLEND 25 . Only 
isolated objects, child objects, and objects on which the deblender gave up are used in constructing 
the galaxy catalog. 



4.2. Star-galaxy classification 

Photo fits each object in all five colors to three model profiles convolved with the locally deter- 
mined PSF, i.e., a point source, an exponential galaxy profile of arbitrary axial ratio, orientation 
and scale length, and a de Vaucouleurs galaxy profile of arbitrary axial ratio, orientation and scale 
length. For each object, the best-fit model in the r' band is noted. The corresponding magnitude 
in the different bands, calculated using the parameters of the best fit model in the r' band is re- 
ferred to as the model magnitude in that band. In principle, it is possible to use the likelihoods of 
these model fits to classify objects into stars and galaxies. We find empirically, however, that the 
star-galaxy classification is better using the following technique. Every object is classified as a star 
or a galaxy in each band, based on the difference between the model and PSF magnitudes 26 . An 
object is classified as a star in any band if the model magnitude and the PSF magnitude differ by 
no more than 0.145 mag, corresponding to the fluxes in the model fit and the PSF fit to the object 
satisfying the relation 

Flux in PSF 1 

Flux in model < 0.875' ^ ' 



Objects touching the edge of a frame will be deblended in the latest versions of Photo. 

26 PSF magnitude is computed for all objects by fitting a two-dimensional PSF model, which is a continuously 
varying function of the position on the frame as determined from bright stars in the data. 
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We have tested the accuracy of this star-galaxy classification algorithm by comparing our 
classifications of objects in SDSS images of the Groth strip against the classifications in the Medium 
Deep Survey catalog (MDS) constructed using WFPC2 parallel images from the Hubble Space 
Telescope (Ratnatunga et al. 1999, and references therein). The SDSS data were taken in May 
2000 (Runs 1468 and 1469). The FWHM of the median seeing was 1.4", slightly worse than the 
FWHM in Runs 752 and 756. We find that our star-galaxy classification correctly reproduces 
the MDS classifications for all objects (60 stars and 68 galaxies) with r* < 20.5. In the range 
20.5 < r* < 21, 51 out of 52 galaxies and 14 out of 14 stars are correctly classified. More details 
and tests of the star-galaxy classification algorithms used in Photo can be found in Lupton et al. 
(unpublished). 



4.3. Bright galaxy sample 

Although the star-galaxy classification algorithm described in the previous subsection yields 
correct classifications in the range 16 < r* < 21 for the purpose of deriving reliable galaxy counts, 
there are reasons to suspect that this algorithm may run into problems at brighter magnitudes. 
For example, saturated pixels and diffraction spikes can compromise the model-fitting algorithm 27 
28 . In constructing the machine galaxy catalog from the photometric catalog, we reject all objects 
that contain saturated pixels (over the entire sky, only a handful of galaxies with active nuclei 
are expected to be saturated in the SDSS images). However, this rejection will also eliminate 
real galaxies that are blended with saturated stars. Hence, in order to calibrate the performance 
of the star-galaxy classification algorithm employed by Photo, and to understand how well the 
deblender measures the magnitudes of galaxies with substantial internal substructure, we have 
visually examined the images of all objects brighter than r* = 16, in the stripe defined by Runs 
752 and 756, covering an area of about 230 deg 2 . At these bright magnitudes and large galaxy 
sizes, we expect our visual classification to be essentially perfect. Hence, we can use the visual 
galaxy catalog to quantify the completeness and the contamination fraction of the galaxy catalog 
constructed by Photo, which we refer to as the machine galaxy catalogue. In this region, Photo 
classifies 27,304 objects as non-stellar (these include saturated stars), and 86,137 objects as stars. 

Table 1 presents the results of comparing the visual and the machine galaxy catalogs at r* < 16. 
Column (1) lists the range of observed magnitudes (not yet corrected for reddening) in the r' band. 
Columns (2) and (3) show the number of objects in that magnitude bin that are classified as galaxies 
by visual inspection and by the machine galaxy catalog, respectively. Column (4) is the number of 
visually selected galaxies that are not included in the machine galaxy catalog. Column (5) shows 
the opposite case; the objects selected as galaxies in the machine galaxy catalog that are not real 
galaxies upon visual inspection. 



Obviously saturated pixels are rejected from fitting, but marginallly saturated pixels may affect the model fitting. 
We're still working on subtracting the diffraction spikes from saturated stars in the current version of Photo. 
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There are 93 objects (5% of the sample) that are classified as galaxies by visual inspection 
but are not present in the machine galaxy catalog (the galaxies in column 4 of Table 1). These 
galaxies are not included in the machine galaxy catalog for the following reasons: (a) 48 galaxies 
are not selected because they are blended with saturated stars, (b) 45 galaxies are not classified 
in the r' band since they contain an interpolated pixel (i.e., a bad column or a cosmic ray) at or 
near the center, a condition which causes Photo to classify an object as UNKNOWN 29 . Among these 45 
galaxies, 43 are correctly classified as galaxies in both g' and i' bands. Hence, these galaxies would 
have been included in the machine galaxy catalog if we had selected as galaxies those objects that 
Photo classifies as galaxies in at least two of the three bands, g', r', and i'. 

Similarly, there are 215 objects that are classified as galaxies in the machine catalog, but are 
not confirmed to be real galaxies upon visual inspection (the objects in column 5 of Table 1). These 
objects are included in the machine galaxy catalog for the following reasons: (a) 108 objects are 
classified as galaxies in the r' band, but are classified as stars in all the other color bands. These 
objects are either single stars with slightly elongated shapes in the r' band, or spurious objects 
caused by satellite trails. None of these objects are real galaxies. They would be excluded from the 
machine galaxy catalog if we required that the object be classified as a galaxy in at least two of the 
three bands g' , r', and i'. (b) 44 objects are double stars misclassified as galaxies, (c) 44 objects 
are redundant children arising from excessive deblending of other bright galaxies already present 
in the visual galaxy catalog, (d) 3 objects are groups of galaxies, whose individual members would 
be fainter than r' = 16 if they had been properly deblended. (e) 4 objects are selected as galaxies 
because of satellite trails and classified as galaxies in at least two of the three colors g', r' and i! . 
(f) 12 objects are ghost images due to reflections of bright stars inside the camera. 

We have re-examined the images of all the 44 spurious children which resulted from excessive 
deblending, since such deblending problems will lead to underestimates of the flux of the parent 
galaxy. We have combined these 44 objects back to 31 parents, and remeasured their Petrosian 
fluxes. We find that the error in magnitude resulting from this excessive deblending ranges from 0.1 
to 1 mag, with a typical value of 0.5 mag. We use the corrected magnitudes for these 31 galaxies in 
the visual galaxy catalog. We note that the frequency of excessively deblending objects decreases 
at fainter magnitudes; in particular, there are no such objects in the 15 < r* < 16 bin, implying 
that excessive deblending is unlikely to be a problem fainter than 16 mag. 

From Table 1, we see that the net statistical completeness of the machine galaxy catalog, 
defined as (Galaxies common to both visual and machine catalogs) /(Galaxies in visual catalog), 
is 95% for the r* < 16 mag sample. However, if we require that all objects in the machine 
galaxy catalog be classified as galaxies in at least two out of the three bands, the completeness 
increases to 97%. The major reason for missing real galaxies (incompleteness) in the machine 
galaxy catalog is that they are blended with saturated stars, while most of the spurious galaxy 
detections (contaminants) are either double stars or shredded pieces of brighter galaxies. The 



Objects with interpolated pixels at the center is properly classified in later versions of Photo. 
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galaxy number counts obtained from the machine galaxy catalog matches the true galaxy counts 
to within 3% for the bright sample with r* < 16. 

The comparison of the Photo outputs to the results of the visual inspection has uncovered a 
number of subtle problems with the photometry processing with Photo version 5.0.2 as of early 
2000. Many of these problems are being fixed for a future version of the pipeline, which should 
have improved performance. 

4.4. Faint Galaxy Sample 

In the magnitude range 16 < r* < 21, we derive the galaxy counts using the machine galaxy 
catalog. At these magnitudes, we construct the machine galaxy catalog from the photometric 
catalog using the following criteria: (a) Only the PRIMARY detection of an object is included in the 
galaxy sample, (b) Based on the experience gained at bright magnitudes, we classify an object as 
a galaxy only if Photo classifies it as a galaxy in at least two of the three photometric bands, g', r', 
and i'. (c) Objects that are blended with saturated stars are not included in the galaxy sample. 

Figure 2 shows the distribution in the g* — r* r* — i* plane of all objects classified as stars 
and galaxies, in fields 251 to 300 of the third CCD column in Run 756 (seeing FWHM is about 
l") 30 . The colors are computed using the model magnitudes of objects in each band. The g* — r* 
color is roughly constant for stars with r* — i* > 0.5 (corresponding to stars cooler than 3500 K), 
as the g 1 and r' bands are dominated by absorption from molecular bands, especially from TiO. 
Hence, the objects with g* — r* > 1.6 that are classified as stars (objects that lie to the right of 
the dashed line in panels (b) and (d)) are quasars, stars with incorrect photometry, or real galaxies 
(including compact galaxies) that are misclassified as stars. Even if we assume that all objects 
with g* — r* > 1.6(1.7) are really galaxies, only 1.4% (1.0%) of galaxies are being misclassified as 
stars in the magnitude range 16 < r* < 20, and 4% (2%) in the range 20 < r* < 21. Similarly, our 
tests of the star-galaxy classification algorithm using the images from the Groth strip (described 
in §4.2) show that at most a few percent of stars are being misclassified as galaxies at r* < 21. We 
have also compared the classifications of objects that lie in the overlap region between Runs 752 
and 756. We found that the two classifications in these runs with different seeing were identical 
for 100% of objects at r* < 20, for 99% of objects in 20 < r* < 20.5, and for 97% of the objects 
in 20.5 < r* < 21. All these tests show that misclassification of objects does not pose a serious 
problem in deriving the galaxy counts at all magnitudes brighter than r* = 21. 



We have not separated quasars from stars. So we expect about 50-100 quasars being included in the star sample. 
We have not studied in this paper the problem of blue compact galaxies. 
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5. Galaxy number counts in the SDSS bands 



Figure 3 shows the variation from one field to another of the number of stars and galaxies, 
averaged over all six CCD columns. The different panels show this variation in different bins in r' 
magnitude. Figure 3(a) shows the objects in Run 756, while Figure 3(b) shows those in Run 752. 
The abscissa in each panel is the field number. Note that a field corresponds to a 10' x 13' region 
of sky, and the data shown are an average over 6 fields across the stripe. The global variation of 
the number of stars (on the left panels) reflects the change in Galactic latitude during the TDI 
scan over 7-8 hours. The top left panel in each figure shows the variation of Galactic latitude 
over the duration of the scan. The number of galaxies varies from field to field due to both small- 
number statistics and small-scale clustering. However, there are no large-scale gradients in the 
galaxy distribution with field numbers in bins brighter than r* = 21, confirming that there is no 
serious contamination of stars to the galaxy counts in the magnitude range of interest. The galaxy 
counts in the faintest magnitude bin at r* > 21 could be systematically affected by galaxies being 
misclassified as stars. In particular, note the drop in the galaxy counts for 21 < r* < 21.5 in Run 
752 for small field numbers, where the seeing is poor, indicating that some fraction may not be 
detected and more importantly some galaxies with low signal to noise ratio could be misclassified as 
stars. At the magnitude where we present the galaxy counts there is no contamination from stars to 
galaxies, despite the fact that there is a higher ratio of stars to galaxies at brighter magnitudes. The 
solid curves in the left panels show the expected variation of star counts from the Galaxy model 
of Bahcall & Soneira (1980) using the color transformations derived from our work on primary 
standard stars (Smith et al. 2001) 



We refer to Chen et al. (2000) and Finlator et al. (2000) for more detailed analysis of star counts. 
These empirical color transformations differ from the synthetic relations derived in Fukugita et al. 
(1996) by less than 0.03 mag in the constant terms. The model predictions and the observed star 
counts agree remarkably well at bright magnitudes (r* < 18.5) for all fields of the TDI scan. For 
18.5 < r* < 20.5, the deviations of the observed star counts from the model predictions become 
increasingly apparent at large field numbers, which are pointing toward low Galactic latitudes in 
the direction of the Galactic center. 

Figure 4 shows the galaxy number counts in 0.5 mag bins in the r' band in the Northern 
equatorial stripe. We have coadded the counts in Runs 752 and 756 for bins at r* < 18. At 
fainter magnitudes, we present the results from Run 756 alone. We are certainly not limited by 
small-number statistics at r* > 18 mag, and therefore we do not include the galaxy counts from 
Run 752, where the seeing is worse. The solid points show the counts from the visual galaxy 
catalog, and the open points show the counts from the machine galaxy catalog. Star counts arc 
also shown as crosses: the sharp drop seen at bright magnitudes (r* < 14.5) is due to the fact 



* 



V + 0.530(5 - V) - 0.075, 

V - 0.438(5 - V) + 0.113. 



(5) 
(6) 
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that we cannot measure stars that saturate in the photometric CCD. Thus, our star counts arc 
reliable in the magnitude range 15 < r* < 21. We deredden the Petrosian magnitudes of galaxies 
using the Galactic extinction map of Schlegel et al. (1998), assuming a ratio of total to selective 
extinction in the r' band R r > = 2.75 31 computed from a standard extinction law (Schlegel et al. 
1998). For most of the sky region we have scanned, extinction is smaller than 0.1 mag (Figure 5 
shows selective extinction E(B — V) along the Celestial Equator). A smooth match between the 
galaxy counts from the visual and the machine galaxy catalogs is apparent; the two samples give 
almost identical galaxy counts in the range of magnitudes 14 < r* < 15. 5. 32 Hence, we derive the 
galaxy counts from the visual galaxy catalog at r* < 15.5, and from the machine galaxy catalog 
at fainter magnitudes. Table 2 presents the galaxy counts in the five SDSS bands u', g', r', i', 
and z' . We list the actual number of galaxies M detected in our sample in each 0.5 mag bin (after 
dereddening the magnitudes of galaxies), together with the effective area in which we searched for 
galaxies in that bin. We note that the effective areas in the magnitude bins from the machine galaxy 
catalog are corrected for regions which are masked to avoid bright stars (amounting to about 0.4% 
of the total area). 

The error bars on the galaxy counts data in Figure 4 include contributions from both Poisson 
noise and large scale structure in the galaxy distribution. The covariance between the galaxy counts 
in two magnitude bins is (Peebles 1993) 



where iVj, Nj are the number of galaxies in bins i and j, cjj, Uj are the areas of the sky from which 
the galaxy numbers are derived, 5ft- is the Kronecker delta symbol, u>ij(9) is the angular correlation 
function between the galaxies in bins i and j, and the double integral extends over every pair of 
points in areas uj-i and LOj that are separated by any angle 9ij . The first term is the contribution from 
Poisson noise and is included only in computing the variance in counts in a magnitude bin. The 
second term is the contribution from large-scale structure in the galaxy distribution. We use the 
covariance matrix of the galaxy counts in different magntiude bins in fitting for the normalization 
of the counts-magnitude relation. 

In any two magnitude bins, we estimate the contribution from the second term as follows. We 
first assume that w(9) can be represented as a power-law of the form w{9) = A(m)9 1 ~ J , where 
7 = 1.8 is the slope of the two point correlation function (Peebles 1993). We describe the details 
of computing the angular correlation function Wij(9) between galaxies in different magnitude bins 



31 Note that Petrosian magnitudes are easy to correct for extinction, unlike, for example, isophotal magnitudes (see 
the discussion in Santiago et al. (1996)). 

32 The visually selected galaxy catalog is no longer complete in the magnitude bin 15.5 < r* < 16 after reddening 
corrections. 
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in Appendix A. We perform the double integral 
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du)idu)j 



(8) 



in a patch of sky 90° x 2.5°, roughly corresponding to the size and geometry of the stripe defined 
by Runs 752 and 756, and use the same value of I\2 for all possible pairs of geometries of the survey 
areas in different magnitude bins. We fix the amplitude of w(9) assuming that w(6 = 0.1°) = 0.1 
for galaxies in the magnitude range 18 < r* < 18.5, approximately the value measured from SDSS 
commissioning data itself (Scranton et al. 2000). 

The line segment in Figure 4 shows a fit to the galaxy counts-magnitude relation expected 
in a homogeneous universe assuming Euclidean geometry for 3 dimensional space (we refer this as 
"Euclidean" geometry in the following text), 



The observed galaxy counts are quite consistent with this line with A r * = 5.99±0.52 (0.5mag)~ 1 deg -2 
(x 2 = 13.9 for 9 degrees of freedom) for 12 < r* < 17 and even beyond this range. For r* < 13, 
the number of galaxies is quite small (N < 20 in each bin), and there are large errors due to both 
Poisson noise and large-scale clustering. 

In order to study the number counts in detail, we plot in Figure 6 the quantity iV(r*)10~ a6<r * _16 ), 
the observed galaxy number counts normalized by the shape of the counts-magnitude relation (here- 
after, called the growth rate) expected for a homogeneous galaxy distribution in a universe with 
"Euclidean" geometry. We also show the predictions of a no-evolution model, taking into account 
both cosmological corrections and K-corrections. The model is based on studies of the statistics of 
galaxy parameters in the B band. We use a mix of E, SO, Sab, and Scd galaxies, in the ratio of 
morphological types given in Fukugita et al. (1998) and use the color shifts and K-corrections given 
in Fukugita et al. (1995). We assume that the luminosity function is characterized by a Schechter 
function with parameters a = —1.1, and Mq = —19.5 + 5 log ft,. The shape of the predicted curves 
depends only weakly on the model parameters. We show the predictions of the no-evolution model 
in three different cosmologies: (fl m ,Q\) = (1,0), (0.3,0), and (0.3,0.7). We also show the pre- 
dictions in a universe assuming Euclidean geometry for three space, but including K-corrections. 
Model curves are normalized to the amplitude A r * derived in equation (10) below. 

The observed galaxy counts are consistent with the predictions of the no-evolution model for 
r* > 16, although the data permit some evolution. At r* < 16, the data mostly lie below the 
predictions of the no-evolution model, but they are mostly within the fluctuations expected from 
Poisson noise and large-scale structure in the galaxy distribution. We will show in § 6 that this 
deficit of galaxies at bright magnitudes is much smaller in amplitude and much more local in 
volume than the local underdensity claimed by the APM group (Maddox et al. 1990). This small 
underdensity toward the Northern Galactic cap is also not seen in our data toward the Southern 
Galactic cap (analyzed in § 8 below). It is therefore likely to be ascribed to a local fluctuation due 
to large scale structure in the galaxy distribution. 



N(m x ) = A A 10 a6(m ^ 16) . 



(9) 
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We repeat the same analysis of the galaxy counts in the u', g', i' and z' passbands. We 
assume the ratio of total to selective extinction to be R' u = 5.16, R' g = 3.79, R[ = 2.09 and 
R' z = 1.48, respectively, taken from a standard extinction law (Schlegel et al. 1998). The results 
are presented in Figure 7 and Table 2. The error bars include contributions from both Poisson 
noise and large scale structure. In each band we estimate w(9) of galaxies in a magnitude bin 
from the corresponding value of w(9) in the r' band, assuming a mean color for galaxies in that 
magnitude bin. The galaxy counts in the g' , i! and z' passbands closely follow the I0 a6m law at 
bright magnitudes (m < 17 — 18), as we have seen with the counts in the r* band. The counts in 
the v! passband show some wiggle for u* < 18 mag, and the local deficit is more conspicuous for 
brighter magnitudes. The reason is not clear to us. 

As is evident from Figure 6, cosmological corrections and K-corrections have a non-negligible 
effect on the galaxy counts at magnitudes even brighter than r' = 16 mag (as seen from the 
difference in shape between the curves corresponding to different cosmologies and the "Euclidean" 
geometry) , although in the data, this effect is masked by fluctuations arising from Poisson noise and 
large-scale structure. To obtain the cosmologically meaningful normalization of the galaxy counts, 
we must take account of the cosmological and K-corrections. The fit must cover a magnitude range 
faint enough that local large scale structure effects are avoided, but not so faint that evolutionary 
effects become important. In Figure 8 (c) we display with solid points the galaxy count data in 
the r' band divided by the no-evolution model prediction (i.e. the "Euclidean" growth factor and 
the cosmological and K-corrections) so that data should fall on a constant line if they are not 
affected by large scale structure and evolution. Here we assume (fi m , ^a) = (0.3, 0.7), although the 
results are insensitive to cosmology (The open points are the counts data from Southern Celestial 
Equator stripe, which are discussed in §8). We fit the data (only those denoted by solid points) 
taking account of Poisson noise and large-scale structure using the fitting range chosen to cover 
the plateau. We have carried out both full correlated and uncorrelated fits. The result of full 
correlated fit, however, turns out to give a rather poor fit to the number count data 33 . Hence, we 
adopt the result from the uncorrelated fit. For example for the r' band, we obtain the coefficient 
of the "Euclidean" growth rate in the bright limit, which we denote by A r *, 

A r * = 11.30 ± 0.75 (0.5 mag)" 1 deg~ 2 . (10) 

where the error includes Poisson noise, large scale structure and the variation induced by the change 
of the faint end of the fitting range by ±1 mag. Figure 8 presents the analysis repeated for all 
passbands (including B band which we discuss in the next section). The la error band is indicated 



33 This implies that the smooth no-evolution model used here is not quite suitable to describe the data (Correlated 
fits are sensitive to the adequacy of the fitting function) . The values of A from the full correlated fit are systematically 
lower by 0.6 — 1.3cr (except for u') than the values from the uncorrelated fit, which we have given in the figure and 
table. The fitting function for v! is apparently inadequate, and we find a larger discrepancy between the two fits. For 
completeness, the values of A from full correlated fit for u* , g* , r* , i* , z* , and B band are 0.91 ± 0.10, 4.73 ± 0.24, 
10.47 ± 0.57, 17.48 ± 0.81, 28.71 ± 1.17, and 2.78 ± 0.19, respectively. 
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with dotted lines in the figure. We obtain the coefficient of the "Euclidean" growth rate in the 
bright limit A in each band, which is given in Table 3, together with the range of magnitudes over 
which it was fit. We could not find a plateau for the z' band (to a lesser extent for the i! band), 
and hence the real error for A(z*) would be larger than is quoted in the table. 

We do not discuss the detail of star counts in this paper. Extensive analyses of the star counts 
in the SDSS bands, and their implications for models of Galactic structure are presented in Chen 
et al. (2000) and Finlator et al. (2000). 

6. Galaxy number counts in the B and I^ia passbands 

The majority of results on galaxy number counts that exist in the literature have used the B (or 
Bj) band. Therefore, in this section, we convert our galaxy counts results to the Johnson-Morgan 
B band photometric system, and determine the normalization of the counts-magnitude relation. 
We adopt the color transformation 

B-g* = 0.482( 5 * - r*) + 0.169, (11) 

derived from preliminary standard star work (Smith et al. 2001) for the Landolt stars (Landolt 
1992) with the standard deviation of 0.03 mag. We find that equation (11) provides a better 
match to the observed B band photometry than an equivalent expression using the u* — g* color 34 . 
We have checked the accuracy of this transformation using colors synthesized from the spectra of 
galaxies. We use the spectroscopic energy distributions of galaxies from the spectrophotometric 
atlas of Kennicutt (1992), and carry out synthetic calculations of u', g', r' and B magnitudes 
(see Fukugita et al. (1995) for more details). The B magnitude calculated via equation (11) from 
synthetic g 1 and r' magnitudes agrees with the synthetic B with an offset of —0.01 to +0.04 mag 
for elliptical and Sb galaxies at zero redshift (+ sign means that the color transformation makes B 
magnitude brighter than the synthetic value). The offset increases to +0.1 mag for late type spirals. 
As we will see in § 9 below, the mean color of our bright galaxy sample is g* — r* ~ 0.7, for which 
the offset nearly vanishes. The color transformation will in general also depend on the redshift of 
the galaxy, due to K-corrections. The accuracy degrades when the 4000A spectral break moves 
deep into the g' band, but we find that the color transformation given by equation (11) works as 
well at z ~ 0.2 as at z ~ 0; for a galaxy at z ss 0.4 the offset increases to about 0.2 mag. We apply 
the color transformation in equation (11) on a galaxy by galaxy basis, and show the results only 
for B < 19.5, as the fainter magnitude bins include a substantial fraction of galaxies with z > 0.3. 

Figure 9 shows the galaxy counts in the transformed B band (after correcting for Galactic 
reddening), and Figure 10 shows the same quantity after normalizing by the "Euclidean" growth 



34 One may consider the transformation to B magnitude via u* and g* , since the B band is located between the u* 
and g* bands. However, this gives a larger offset except for late type spiral galaxies; for example, for g — r' ~ 0.7 
the offset is 0.06 mag. This offset results from the difference of u' band spectral features of stars and galaxies. 
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factor of 10 6m . The latter fi gure also shows the galaxy counts data in the B band from the APM 
survey (Maddox et al. 1990), and from the surveys of Gardner et al. (1996) and Bertin & Dennefeld 
(1997), together with the predictions of the no-evolution model. The photographic bj band data of 
Maddox et al. (1990) and Bertin & Dennefeld (1997) have been converted to B band by Gardner 
et al. (1996). Note that each of these studies uses a different method to determine the magnitudes 
of galaxies, although there are attempts to make corrections to obtain pseudo-total magnitudes; 
we refer the reader to the respective papers for details. 

The shape of our galaxy counts fainter than B = 16.5 agrees well with those of Gardner et 
al. (1996) and Bertin & Dennefeld (1997), but our normalization is about 15% higher than theirs. 
On the other hand, the APM data lie substantially below our counts brighter than 18.5 mag. In 
particular, in the magnitude range 15 < B < 18, our galaxy counts are on average a factor of 
two larger than the counts inferred from the APM survey. The shape of our counts-magnitude 
relation is also consistent with the predictions of the no-evolution model in the magnitude range 
17 < B < 19, while some deficit is seen in magnitudes brighter than B = 16 arising from local 
large scale structure in the galaxy distribution. Since the depth corresponding to the deficit is only 
150/i _1 Mpc^" 1 and our survey is only in a small patch on the sky, it appears that the deficit brighter 
than B = 16 is consistent with the amplitude of large-scale structure fluctuations. We need not 
invoke an overall local underdensity of galaxies over a large solid angle. This is corroborated by 
the analysis for the Southern stripe, as we discuss later. We note that the APM data remain flat 
over the range from 15.5 to 19 mag, while cosmological effects and K-correction should make the 
counts decrease significantly over this magnitude range; Maddox et al. (1990) interpret this as due 
to evolution. Our data do not show evidence for such a rapid 'evolution' of the galaxy population 
to z = 0.2 (approximately the redshift of an L* galaxy at B = 19). 

The star counts plotted in Figure 9 agree well with the predictions of the Bahcall-Soneira 
model in the range 15 < B < 19. The agreement is a qualitative confirmation that our color 
transformation in equation (11) works well for stars. 

We have also calculated the galaxy number counts in the isi4 passband (F814W) — the band 
in which much work has been done with the Hubble Space Telescope. From our standard star work 
(Smith et al. 2001), we find for the transformation to the Cousins' I c band that 

I c - i* = -0.205(r* - i*) - 0.382, (12) 

with the standard deviation of 0.02 mag. We further apply color transformation from I c to I814: 
^814 = Ic + 0.04 according to Fukugita et al. (1995). This constant depends little on galaxy colors. 

Figure 11 presents our galaxy counts in the isi4 band for igi4 < 20 mag. We also plot the 
data from Gardner et al. (1996) and Postman et al. (1998). The three independent data sets show 
good agreement up to an offset of about 0.05 mag, which might be ascribed to systematic errors in 
the different bands or to different definitions of the magnitudes of objects. Our numerical data for 
the galaxy counts are given in Table 4 for both B and ^814 bands. 
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7. Local luminosity density of the universe 

The normalization of the number counts A of eq. (9) is related to the parameters of the 
Schechter luminosity function (Schechter 1976) as 

.4 = 0.7046^ (d* 6 ) 3 r(a + 5/2)0*, (13) 

where d* = lOpc x io a2 ( 16 - M *), and lo is the solid-angle coverage of the sample. When we take ac- 
count of cosmological effects and K-corrections, A must be replaced with the coefficient A obtained 
after corrections as we saw in §5. We then compute 0* using equation (13). 

Let us calculate the luminosity density of the universe in the B band, which has long suffered 
from significant uncertainties. Since we do not determine M* and a in the B band from our own 
data, we use the measurements from the literature (Efstathiou et al. 1988; Loveday et al. 1992; Zucca 
1997; Ratcliffe 1998; Folkes et al. 1999). These parameters lie in the range -19.7 < M* B < -19.4 
and —1.0 > a > —1.2 35 , where Bj = B + OA is applied to convert into the B band when necessary. 
From A B = 2.98 ± 0.17 (B = 12 - 19.5), we obtain (f>* = 2.05 ± 0.12+°;^. Here the first error 
stands for that of A and the second arises from uncertainties in the Schechter parameters. 

We then calculate the luminosity density of the universe in the B band. Note that the fractional 
error in luminosity density (£) is much smaller than that in </>*, because the luminosity density 
goes as £ ~ (p*L*T(a + 2) while the number counts go roughly as N(m) ~ <f>* L* 3 / 2 T(a + 5/2), so 
that errors in L*, a and 4>* largely cancel. We find that the error in C produced by uncertainties in 
the Schechter parameters is reduced to 3%. The error in C is dominated by the error in the global 
number counts. We obtain 

C B = (2.41 ± 0.39) x 10 8 L o /i(Mpc)- 3 . (14) 

as our best estimate. Here we have applied a 3% upward shift for the expected offset between our 
Petrosian flux and the total flux (see the discussion in § 3), and included a 5% error from provisional 
photometric calibration and uncertainties in the color transformation equations. 

We also calculate the luminosity density of the universe in the five SDSS passbands u' , g', r', i' 
and z' , in a similar manner, adopting the values of the Schechter luminosity function parameters 
derived from the SDSS commissioning data itself (Blanton et al. 2000). Table 5 presents our 
estimates of cf)* x and C\ in each band (in units of the solar luminosity in the respective band). 
We adopt M u #(©) = 6.38 M g ,(&) = 5.06, M r ,(&) = 4.64, M V {Q) = 4.53 and M Z ,(Q) = 4.52 
which roughly correspond to (B — V)q = 0.63 (Fukugita et al. unpublished). Here we assume the 
covariance of the errors of the Schechter fit, A(Mb ■ aj/^AMe) 2 • (Aa) 2 w 1. The values of <f>* x 



35 Folkes et al. (1999) give a significantly steeper than this range. Their Schechter fit, however, poorly represents 
the step-wise maximally likelihood result, and the faint end slope is significantly flatter than their a indicates. 
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and C\ are consistent with those of Blanton et al. (2000) within our la error for r',i'andz'. Our 
values are higher for u'andg', but within 2a. 

8. Galaxy number counts in the Southern equatorial stripe 

We have also analyzed the galaxy counts in an additional equatorial stripe toward the Southern 
Galactic cap. These data were taken on September 19 and 25, 1999 (SDSS Runs 94 and 125). The 
first run scanned the sky for 5.3 hours from a = 336.7° to 56.28°, with declination extending from 
— 0.°93 to +1.°17. The second run scanned for 5.7 hours from a = 350. °6 to 76. °7 with declination 
varying from — 1.°14 to +0.°96; the two runs together make a filled stripe. The median seeing 
FWHM in the r' band was 1.4" in Run 94 and 1.7" in Run 125, but the telescope was not well 
collimated during these runs. The total area covered in this Southern equatorial stripe is 210 deg 2 . 
We do not construct the visual galaxy catalog in this region, but instead use the machine galaxy 
catalog itself over the entire range of magnitudes 11 < r* < 21. For r* > 18, we present the galaxy 
counts from Run 94 alone, since it has data taken under better seeing. We present the galaxy 
counts in the u', g', r', i' and z' bands in Figure 13 and Table 6. 

Figure 12 shows the galaxy number counts in r' band normalized by the "Euclidean" growth 
factor of 10°- 6m , in both the Northern (solid points) and the Southern (open points) equatorial 
stripes. The error bars on the galaxy counts in both stripes include contributions from both 
Poisson noise and large scale structure. The number counts from these two essentially independent 
regions of the sky are identical for r* > 16 mag. At brighter magnitudes, the counts in the Southern 
stripe are systematically larger than those in the Northern stripe. The Southern stripe crosses the 
southern extension of the Pisces-Perseus supercluster; the enhanced counts are due to this well- 
known overdensity. The curves represent the predictions of the no-evolution model we described 
earlier. 

The degree of disagreement in the galaxy counts at r* < 16 between the Northern and the 
Southern stripes indicates that the deficit of galaxies seen in the Northern stripe is a local effect 
that depends on direction in the sky; i.e., the local universe to 100—200 Mpc is clumpy. In the 
magnitude range 14 < r* < 15.5, there are N\ = 1010 galaxies in the Northern stripe, in an area 
of uoi = 230 deg 2 . The rms fluctuation in this number due to large scale structure is computed 
using equation (7) to be a = 170. In the same magnitude range, there are N2 = 1259 galaxies 
in the Southern stripe, in an area of u)2 = 210 deg 2 . The difference in counts between the two 
stripes is therefore \N\ — iV^i/u^l = 369 — plausibly a 2a fluctuation due to large scale structure 
in the galaxy distribution. Thus, the difference in bright galaxy counts between the Northern and 
Southern equatorial stripes can be explained as arising from a large scale structure fluctuation in 
the galaxy distribution. This supports our earlier observation in § 6 that the deficit of galaxies 
at bright magnitudes in the Northern stripe is due to large scale structure fluctuations. We have 
to wait, however, for the end of the survey to draw definitive and quantitative conclusions on this 
point, as we quantify large scale structure over a much larger sample. 
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We refer to Figure 8 above for a comparison of the counts in the Northern and Southern 
equatorial stripes for the other color bands. 

9. Color distributions of galaxies 

Figures 14 and 15 show the distribution of g* — r* and r* — i* colors of all galaxies in the third 
CCD column of Run 756 as a representative sample. Dereddened Petrosian magnitudes are used 
to calculate colors, and color distributions are shown according to magnitude bin in r*. Note that 
the Petrosian magnitude in each color band is computed using the Petrosian aperture calculated 
in the r' band. Therefore, the colors are measured through a consistent aperture in all bands. The 
mean (/x) and rms scatter (a) are indicated in each panel. 

The mean of g* — r* colors for the brightest bin is 0.71, which is slightly (0.04 mag) redder 
than expected from a spectroscopic synthesis calculation for a standard morphology mix (Fukugita 
et al. 1995). Figure 14 shows that the mean g* — r* colors of galaxies become systematically bluer 
at bright magnitudes; the mean value of g* — r* gradually shifts to 1.03 at 20 < r* < 21. The mean 
r* — i* color of galaxies (0.36 in Figure 15) is consistent with what is expected from a synthesis 
calculation. This color stays roughly constant over the entire range of magnitudes 15 < r* < 21, 
as expected from the spectroscopic energy distributions of normal galaxies. 

The dispersions in both the color distributions increase towards faint magnitudes. The scatter 
in color in bin 20 < r* < 21 is 0.52 mag in g* — r*, and 0.35 mag in r* — i*. These scatters are much 
larger than those expected from photometric errors alone; the photometric error in colors in this 
magnitude range is about 0.10 mag in g* — r*, and 0.07 mag in r* — i*, both as quoted by Photo 
and from tests using simulated images. This broadening of the color distributions is expected from 
the K-corrections, which become increasingly different between galaxies of different morphological 
types with increasing redshift, and also from the broader range of redshifts that contribute to the 
colors in a given magnitude bin at fainter magnitudes. 

10. Conclusions 

We have presented the number counts of galaxies in five SDSS bands, u', g', r', i' and z', in 
the magnitude range 12 < r* < 21. We have first examined in detail the photometric catalog 
derived from automated software, Photo. We have carefully analyzed the bright galaxy sample by 
comparing the galaxy catalog produced by Photo against that produced by eye inspection. We 
have shown that the statistical completeness of the machine generated galaxy catalog is 97% for 
the sample with r* < 16 mag. However, 5% of objects in the machine galaxy catalog up to r* = 16 
mag are not real galaxies. Toward the faint end, the imaging data are complete to r* = 22.5 mag, 
and our star-galaxy classification is 100% accurate to r* = 20 mag. The misclassification rate is 
still a few percent at r* = 21 mag. We expect that the performance will be improved in the future 
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version of Photo. 

We have then presented the number counts of galaxies in five SDSS bands, u', g', r', i' and z', in 
the magnitude range 12 < r* < 21. The galaxy number counts are derived from two independent 
stripes of imaging data along the Celestial Equator, one each toward the North and the South 
Galactic cap, each covering about 200 deg 2 of sky. The galaxy counts from both stripes agree very 
well in the magnitude range 16 < r* < 21, implying that we are sampling galaxies in a volume that 
is large enough to be a fair sample of the universe at these magnitudes. At brighter magnitudes 
(14 < r* < 16), the galaxy counts from the two stripes differ by about 30%. This difference is 
consistent (at the 2a level) with the fluctuations arising from large scale structure in the galaxy 
distribution. Additionally, we have used empirically determined color transformations to translate 
galaxy counts in the SDSS color bands into those in the B and I^ia bands. We obtained the B 
band local luminosity density, Cb = 2.4 ± 0.4 x 10 8 L Q /iMpc -3 , for a reasonably wide range of 
parameters of the Schechter luminosity function in the B band. 

At the bright magnitudes where cosmological and evolutionary corrections are relatively small, 
the shape of the galaxy number counts-magnitude relation is well characterized by N(m\) oc 
lQ0.6m A ^ foe relation expected for a homogeneous galaxy distribution in a universe with "Euclidean" 
geometry. We have used the amplitude of this relation, A\, to determine the luminosity density of 
the universe at zero redshift in all the five SDSS bands and in the B band (see Table 5). 

The results presented in this paper are based on only a fraction of the imaging data taken 
during the commissioning phase of the SDSS. However, given the high rate of imaging of the 
SDSS (of about 20 deg 2 per hour), we already have large sufficient sky coverage from only a few 
nights of observations to place interesting constraints on such fundamental quantities as the counts- 
magnitude relation and the luminosity density of the universe. The uniform photometric accuracy 
and large dynamic range of the CCDs, coupled with the high rate of imaging of the SDSS will soon 
enable studies of both the properties and the large scale distribution of galaxies at unprecedented 
accuracy. 

The errors in the galaxy number counts and the luminosity density of the universe in different 
bands derived in this work are still dominated by fluctuations due to large scale structure. The 
full SDSS survey will increase the sky coverage by a factor of 25, and more importantly enable the 
study of field to field variation of bright galaxy counts. Since the fractional error in galaxy counts 
due to large scale structure is inversely proportional to sky coverage, we can expect to reduce this 
error (and the error in the luminosity density of the universe) by a similar factor. 
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A. Angular cross-correlation between galaxies in different magnitude bins 

The relation between the angular two-point auto-correlation function w a {9) computed from 
a flux-limited galaxy photometric catalog, and the spatial two-point correlation function £(r), is 
(Limber 1953; Peebles 1993) 

„, ff a \ I J(nr 2 ) 2 rfr 1 (ir 2 ^(r 1 )^(r 2 )e(ri2) , . 

Wa{L e) = [jv^wF ' ( } 

where / is the flux limit of the photometric catalog, r^i = \/r{ + r\ + 2rir 2 cos(0) is the distance 
between the points r\ and r 2 , and the selection function ip(r) is related to the galaxy luminosity 
function <fi{L) by 

poo 

M) = / <f>(L)dL. (A2) 

J 4-Kr 2 f 

Note that in deriving equation (A2) we have assumed that the galaxy clustering is independent of 
luminosity. The angular auto-correlation function w a (f±, / 2 ,6*) of galaxies with fluxes in the range 
/ 2 < / < fi is similar to equation (Al), but with the selection function defined by 



/>47rr 2 /i 

M) = / 4>(L)dL. (A3) 

J47rr 2 /2 



The two-point angular cross-correlation between galaxies in two different magnitude bins cor- 
responding to flux ranges / 2 < / < /i and f± < f < is 

in(f f f f q , _ ! !{nr 2? dndr 2 Mr i)Mr2)^{r 12) . . 

w c {h, hi J3, H, V) r a 1 — r 2a 1 ■ \ A4: ) 

J r^dripi J r z drip2 

The two selection functions are 

rA-nr 2 fi i-A-wr 2 

Mr) = / <KL)dL, Mr) = / <f>(L)dL. (A5) 

1 /4 7rr 2j 2 J47rr 2 /4 

For equally spaced magnitude bins, (/4// 2 ) = (/3//1), 

Mr) =Mr\[h[T2)- (A6) 
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For a power-law correlation function of the form i{r\ 2 ) = (f^r) 7 > equation (A4) becomes 



Wc(fl,f2,f3,U,0) 



r o 1 1 ( r i r 2) 2 dridr 2 ipi(ri)'ipi(r2y/U/h)r 



12 



/ r 2 dripi(r) J r 2 drip 1 (r^/f 4 /f 2 ) 
Under the small-angle approximation <C 1), 

r 12 = y/(r 2 - r{) 2 +r 1 r 2 9 2 . 

The numerator in equation (A7) is dominated by small values of r\ 2 . So, defining 

r 2 - n 



v = 



ro 



(A7) 



(A8) 



(A9) 



and setting n = r 2 = r in the numerator of equation ( A7) , we get 

3/2 r /• i -2 



"', (J\.- k- h- J\-0) i\) ( jr 



r dripi(r) 



J r^r)Mr\[hFk)dr J — 



dv 



\v 2 + r 2 6 2 ] l/2 

(A10) 

where the integral over v ranges from — oo to +oo (since both r\ and r 2 vary between to +oo). 
Now, 



/ 



dx 



(a 2 + x 2 ) 1 ' 2 
Therefore, equation (A10) becomes 



dx 



(1+x 2 ) 



2^7/2 



w c {h, h, h, h,9) = 



H O 1 ^ (^) 3/2 I r5 ~ ldr ^ r )Mrh/k) 
7 \k) [/ r 2 dr^i(r)] 2 



where 



dx 



[J* r 2 dripi(r)Y 
(-l/2)!( 7 /2-3/2)! 

(l +x 2) 7 /2 (7/2-1)! 

Since w a (fi, f 2 ,0) = w c (fi, f 2 , fi, f 2 , 9), equation (A12) can be written as 

„,(f f f f 0)-m(f f 0) (fA 3/2 f r5 -"Mr)MrVf^k)dr 



/+00 
-00 (1 



k 



f r 5 -T [ipi{r)] z dr 



(All) 



(A12) 



(A13) 



(A14) 



We use this relation between w c (fi,f 2 , k, k, 6) and w a (fi, f 2 , 6) to construct the covariance matrix 
of the galaxy counts in different magnitude bins. We evaluate the integrals in equation (A14) using 
the luminosity function in the appropriate SDSS passband derived by Blanton et al. (2000). 
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Fig. 1. — Field to field variation of seeing (FWHM) in the r' band in the third CCD column, during 
TDI observations: (a) Run 752, extending for 7 hours, (b) Run 756, extending for 8 hours. 
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Fig. 2. — Distribution in the g* — r* , r* — i* plane of all objects classified as galaxies and stars in 
fields 251 to 300 of the third CCD column of Run 756. (The star sample contains 50-100 quasars.) 
The inner parts of the distributions are shown as contours, linearly spaced in the density of objects 
in the color-color plane. Panels (a) and (b) show the distribution of objects classified as galaxies and 
stars, respectively, in the magnitude range 16 < r* < 20. Panels (c) and (d) show the corresponding 
objects in the magnitude range 20 < r* < 21. Objects in panels (b) and (d) with g* — r* > 1.6 
(shown by the dashed line) are real galaxies that are misclassified as stars, quasars, or stars with 
incorrect photometry. 
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(a) Run 756 
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(a) Run 756 — continued 
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Fig. 3. — Field to field variation of the number of stars and galaxies in (a) Run 756 and (b) Run 
752, averaged over all six CCD columns. The solid curves in the left panels show the variation of 
star counts predicted by the Bahcall-Soneira model. The top left panel also shows the variation of 
Galactic latitude over the duration of this run. 
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(b) Run 752 — continued 




20.5 < r* < 21 



400 600 800 
Field 



Galaxies 
19.5 < r* < 20 



I 1 1 1 I 1 1 1 I 

20 < r" < 20.5 



20.5 < r* < 21 



200 400 600 800 
Field 



Fig. 3.- 



continued 



-28- 




10 15 20 

r* [mag] 

Fig. 4. — Number counts of galaxies as a function of magnitude in the r' band. Solid points show 
the galaxy counts from the visually inspected sample, and open points from the machine-selected 
sample. The error bars include contributions from both shot-noise and large scale structure (see text 
for details). The line segment shows the counts-magnitude relation expected for a homogeneous 
galaxy distribution in a universe with "Euclidean" geometry: N(r*) = ^4 r *10 a6r *. The crosses 
show the observed star counts (small crosses at r* < 14.5 show the data where stars saturate in 
the SDSS images, and therefore suffer from incompleteness). The solid curve shows the prediction 
of the Bahcall-Soneira model. 
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Fig. 5. — Selective extinction E(B — V) from Schlegel et al. (1998) along the Celestial Equator. 
The solid curves show the ranges included in Runs 752/756 (North) and Runs 94/125 (South). 



-30- 




Fig. 6. — Galaxy counts-magnitude relation in the r* band, normalized by the expected growth 
rate in an "Euclidean" universe i.e. N(r*) x lo[~ a6(r *~ 16 )] . Points show the observed galaxy counts, 
while the curves show the predictions of a no-evolution model in three different cosmologies. All 
the model curves are normalized to the amplitude A r * derived using equation (10). 
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Fig. 7. — Number counts of galaxies as a function of magnitude in the five color bands u', g', r', %' 
and z' for the Northern equatorial stripe. Solid points show the galaxy counts from the visually 
inspected sample, and open points from the machine-selected sample. The error bars include con- 
tributions from both shot-noise and large scale structure (see text for details). The line segment 
shows the counts-magnitude relation expected in a homogeneous universe with "Euclidean" geom- 
etry: N(m) = A m l0°- 6m . The crosses show the observed star counts (small crosses show the data 
where stars saturate in the image, and therefore suffer from incompleteness), and the solid curve 
shows the prediction of the Bahcall-Soneira model. 
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Fig. 7. — continued 
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Fig. 8. — Galaxy counts-magnitude relation in the u', g', r', i' z' , and B bands, normalized by 
the no-evolution model prediction N(r f ) no cvo \. Solid points show the observed galaxy counts from 
the Northern Equatorial stripe, and open points are those from the Southern Equatorial stripe as 
discussed in §8. The solid lines correspond to A\* and dotted lines show one sigma error ranges. 
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Fig. 8. — continued 
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Fig. 9. — Number counts of galaxies as a function of magnitude in the B band, using the color 
transformation given in equation (11). Solid points show the galaxy counts from the visually 
inspected sample, and open points from the machine-selected sample. The error bars include con- 
tributions from both shot-noise and large scale structure (see text for details). The line segment 
shows the counts-magnitude relation expected in a homogeneous universe with "Euclidean" geom- 
etry: N(B) = AbW 0SB . The crosses show the observed star counts (small crosses show the data 
where stars saturate in the image, and therefore suffer from incompleteness), and the solid curve 
shows the prediction of the Bahcall-Soneira model. 
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B [mag] 

Fig. 10. — Galaxy counts-magnitude relation in the B band, normalized by the expected growth 
rate in an "Euclidean" universe i.e. N(B) x 10[~°- 6 ( B ~ 1G ^ . Solid points show the observed galaxy 
counts from this work, while the curves show the predictions of a no-evolution model in three 
different cosmologies. All the model curves are normalized to the amplitude Ab of the data. The 
open circles, crosses and open squares show the data on galaxy counts in B band from three previous 
studies. 
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Fig. 11. — Number counts of galaxies as a function of magnitude in the Jsi4 band, using the color 
transformation given in equation (12) . Solid circles show the galaxy counts from the visually in- 
spected sample, and open circles from the machine-selected sample. The error bars include contribu- 
tions from both shot-noise and large scale structure. The line segment shows the counts-magnitude 
relation expected in a homogeneous universe with "Euclidean" geometry: N(I) = A/10 - 67 . The 
crosses show the observed star counts. The squares and triangles show galaxy counts from previous 
studies. 
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Fig. 12. — Galaxy counts-magnitude relation in the r* band, normalized by the expected growth 
rate in an "Euclidean" universe i.e. N(r*) x 10[~ a6 ( r * _16 )l. Solid points show the galaxy counts 
from the Northern equatorial stripe, and open points show the counts from the Southern equatorial 
stripe. The curves show the predictions of a no-evolution model in three different cosmologies. The 
model curves are normalized to the amplitude A r * of the data from the Northern equatorial stripe. 
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Fig. 13. — Number counts of galaxies as a function of magnitude in the five color bands u', g' , r', i' 
and z 1 for the Southern equatorial stripe. Open points show the galaxy counts from the machine- 
selected sample. The error bars include contributions from both shot-noise and large scale structure 
(see text for details). The line segment shows the counts-magnitude relation expected in a homo- 
geneous universe with "Euclidean" geometry: N(m) = A m W°- 6m . The crosses show the observed 
star counts (small crosses show the data where stars saturate in the image, and therefore suffer 
from incompleteness), and the solid curve shows the prediction of the Bahcall-Soneira model. 
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Fig. 13. — continued 
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Fig. 14. — Distribution of g* — r* colors of all galaxies in the third CCD column of Run 756. The 
colors are computed from the dereddened Petrosian magnitudes of galaxies. Panels (a) to (f) show 
this distribution in different bins of r' Petrosian magnitude. The mean (^) and scatter (<r) of the 
distribution are listed in each panel. 
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Fig. 15. — Distribution of r* — i* colors of all galaxies in the third CCD column of Run 756. The 
colors are computed from the dereddened Petrosian magnitudes of galaxies. Panels (a) to (f ) show 
this distribution in different bins of r* Petrosian magnitude. The mean (//) and scatter (a) of the 
distribution are listed in each panel. 
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Table 1. Comparison of the visual and machine galaxy catalogs at bright magnitudes. 



Magnitude range 


Visual counts 


Machine counts 


Visual only 


Machine only 


9.5-10.0 


1 





1 





10.0-10.5 














10.5-11.0 


1 


1 








11.0-11.5 


4 


1 


3 





11.5-12.0 


4 


5 





1 


12.0-12.5 


4 


3 


1 





12.5-13.0 


10 


14 





4 


13.0-13.5 


30 


32 


3 


5 


13.5-14.0 


58 


75 


4 


21 


14.0-14.5 


101 


114 


12 


25 


14.5-15.0 


257 


279 


15 


37 


15.0-15.5 


517 


548 


24 


55 


15.5-16.0 


987 


1024 


30 


67 


Total 


1974 


2096 


93 


215 



-44- 



Table 2. Number counts of galaxies in the u', g', r', i', and z' color bands for the Northern 
equatorial stripe. N is the number of galaxies in the sample, and N is the normalized counts in 
units of counts per 0.5 mag per deg 2 . Results flagged with a V are based on the visual galaxy 

catalog (§ 4), while the rest are based on the machine galaxy catalog. The galaxy counts at 
magnitudes fainter than 18 are based on Run 756 alone. 



Mag range 
(mag) 


Area 
(deg 2 ) 








N g , 


N r * 


N r * 




JV 4 . 






11.0-11.5 


228v 


1 


0.004 


3 


0.013 


4 


0.018 


3 


0.013 


2 


0.009 


11.5-12.0 


228v 


1 


0.004 


3 


0.013 


3 


0.013 


3 


0.013 


3 


0.013 


12.0-12.5 


228v 


1 


0.004 


4 


0.018 


6 


0.026 


16 


0.070 


26 


0.114 


12.5-13.0 


228v 


3 


0.013 


4 


0.018 


23 


0.101 


28 


0.123 


33 


0.144 


13.0-13.5 


228v 


3 


0.013 


18 


0.079 


28 


0.123 


61 


0.267 


77 


0.337 


13.5-14.0 


228v 


5 


0.022 


28 


0.123 


73 


0.320 


97 


0.425 


145 


0.635 


14.0-14.5 


228v 


11 


0.048 


59 


0.258 


118 


0.517 


244 


1.07 


352 


1.54 


14.5-15.0 


228v 


12 


0.053 


96 


0.420 


315 


1.38 


523 


2.29 


710 


3.11 


15.0-15.5 


228v 


24 


0.105 


233 


1.02 


577 


2.53 


924 


4.05 








283 


















1819 


6.44 


15.5-16.0 


228v 


53 


0.232 


452 


1.98 
















283 










1513 


5.35 


2597 


9.20 


3734 


13.2 


16.0-16.5 


283 


222 


0.784 


1170 


AAA 


3050 


10.8 


4806 


17.0 


7268 


25.7 


16.5-17.0 


283 


555 


1.96 


2290 


8.1 


5607 


19.8 


9337 


33.0 


13601 


48.1 


17.0-17.5 


283 


1393 


Am 


4312 


15.3 


10532 


37.3 


16474 


58.3 


24426 


86.4 


17.5-18.0 


283 


3121 


11.0 


7939 


28.1 


18239 


64.5 


29553 


104.6 


43837 


155.1 


18.0-18.5 


149 


2823 


18.9 


6976 


46.8 


16890 


113.2 


27232 


182.5 


41673 


279.4 


18.5-19.0 


149 


4898 


32.9 


11653 


78.1 


29443 


197.4 


47046 


315.4 


73714 


494.1 


19.0-19.5 


149 


8105 


54.4 


19359 


129.8 


50206 


336.5 


77322 


518.3 


121858 


816.9 


19.5-20.0 


149 


14799 


99.3 


31606 


211.9 


80324 


538.4 


121858 


816.9 


151616 


1016 


20.0-20.5 


149 


28359 


190.3 


51000 


341.8 


123213 


825.9 


166314 


1114 


117358 


786.7 


20.5-21.0 


149 


60618 


406.8 


82614 


553.8 


181428 


1216 


145706 


976.7 


64083 


429.6 


21.0-21.5 


149 


111725 


749.8 


131641 


882.4 


242743 


1627 


65955 


442.1 


31661 


212.2 



-45- 



Table 3. Normalization coefficients of the I0 a6m law, including the cosmological and 

K-corrections 



Band 


Fitting range 




(A) 


(mag) 


[(0.5mag) -1 deg~ 2 ] 


* 

u 


12.0-19.5 


1.52 ±0.18 


9* 


12.0-19.0 


4.95 ± 0.29 


r* 


12.0-18.0 


11.30 ±0.75 


i* 


11.5-18.0 


18.66 ±0.90 


z* 


11.5-18.0 


30.17 ±2.40 


B 


12.0-19.5 


2.98 ±0.17 



-46- 



Table 4. Number counts of galaxies in the B and isi4 bands in the Northern equatorial stripe 

(see Table 2 for explanations) 



Mag range 
(mag) 


Area 
(deg 2 ) 


Mb 


N B 




N hl4 


11.0-11.5 


228v 





0.000 


1 


0.004 


11.5-12.0 


228v 


4 


0.018 


11 


0.048 


12.0-12.5 


228v 


1 


0.004 


27 


0.118 


12.5-13.0 


228v 


6 


0.026 


58 


0.254 


13.0-13.5 


228v 


4 


0.018 


86 


0.377 


13.5-14.0 


228v 


21 


0.092 


208 


0.911 


14.0-14.5 


228v 


24 


0.105 


455 


2.00 


14.5-15.0 


228v 


63 


0.276 


854 


3.74 


15.0-15.5 


228v 


108 


0.473 








283 






2377 


8.41 


15.5-16.0 


228v 


237 


1.04 








283 






4347 


15.5 


16.0-16.5 


283 


608 


2.15 


8501 


30.1 


16.5-17.0 


283 


1230 


4.35 


15236 


53.9 


17.0-17.5 


283 


2423 


8.57 


27316 


96.7 


17.5-18.0 


283 


4526 


16.0 


48717 


172.4 


18.0-18.5 


149 


4167 


27.9 


44432 


297.8 


18.5-19.0 


149 


7192 


48.2 


74565 


499.8 


19.0-19.5 


149 


11569 


77.5 


118000 


791.0 


19.5-20.0 


149 


19017 


127.4 


157416 


1055 



-47- 



Table 5. Luminosity density of the universe in different bands. 



Band 


M* 




ft 

10~ 2 h 3 Mpc -3 


M 


10 s L /iMpc- 3 


u* 


-18.34 ±0.08 


-1.35 ±0.09 


5.67 ± 1.05 


6.38 


6.30 ±0.85 


# 

9 


-20.04 ± 0.04 


-1.26 ±0.05 


1.81 ±0.16 


5.06 


2.53 ±0.20 


r* 


-20.83 ± 0.03 


-1.20 ±0.03 


1.40 ±0.11 


4.64 


2.60 ±0.21 


i* 


-21.26 ±0.04 


-1.25 ±0.04 


1.27 ±0.10 


4.53 


3.31 ±0.23 


z* 
B 


-21.55 ±0.04 
-19 60 +a2 ° 


-1.24 ±0.05 
-1.1 ±0.1 


1.38 ±0.14 

Z - UO -0.30 


4.52 
5.46 


4.60 ± 0.43 

241 +o.46 



-48- 



Table 6. Number counts of galaxies in the u', g', r', i', and z' bands in the Southern equatorial 

stripe (see Table 2 for explanations) 



Mag range 
(mag) 


Area 
(deg 2 ) 












N r , 


K* 






N z * 


14.0-14.5 


209 


6 


0.029 


91 


0.436 


196 


0.939 


352 


1.686 


453 


2.17 


14.5-15.0 


209 


39 


0.187 


149 


0.714 


413 


1.98 


578 


2.77 


734 


3.52 


15.0-15.5 


209 


108 


0.517 


322 


1.54 


650 


3.11 


968 


4.64 


1381 


6.62 


15.5-16.0 


209 


228 


1.09 


576 


2.76 


1157 


5.54 


1842 


8.83 


2731 


13.08 


16.0-16.5 


209 


475 


2.27 


984 


4.71 


2262 


10.8 


3477 


16.65 


5251 


25.1 


16.5-17.0 


209 


843 


4.03 


1811 


8.68 


3962 


18.9 


6461 


30.9 


10063 


48.2 


17.0-17.5 


209 


1733 


8.29 


3274 


15.7 


7525 


36.0 


12245 


58.7 


18648 


89.3 


17.5-18.0 


209 


3197 


15.3 


5501 


26.4 


13823 


66.2 


22653 


108.5 


34586 


165.7 


18.0-18.5 


100 


2056 


20.6 


4502 


44.8 


11724 


116.8 


18947 


188.7 


29486 


293.7 


18.5-19.0 


100 


3426 


34.3 


7786 


77.5 


20435 


203.5 


32309 


321.8 


52662 


524.5 


19.0-19.5 


100 


5752 


57.5 


13698 


136.4 


34895 


347.5 


52893 


526.8 


84069 


837.3 


19.5-20.0 


100 


11178 


111.8 


22262 


221.7 


55831 


556.0 


81855 


815.2 


93168 


927.9 



